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Abstract 

Based on the work of Gorling and that of Levy and Nagy, density-functional formalism for many 
Fermionic excited-states is explored through a careful and rigorous analysis of the excited-state 
density to external potential mapping. It is shown that the knowledge of the ground-state density 
is a must to fix the mapping from an excited-state density to the external potential. This is 
the excited-state counterpart of the Hohenberg-Kohn theorem, where instead of the ground-state 
density the density of the excited-state gives the true many-body wavefunctions of the system. 
Further, the excited-state Kohn-Sham system is defined by comparing it's non-interacting kinetic 
energy with the true kinetic energy. The theory is demonstrated by studying a large number of 
atomic systems. 
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I. INTRODUCTION 



Hohenberg-Kohn-Sham(HKS) density functional theory(DFT)i'2i^ has been most widely 
used to investigate the electronic structures of many-electron systems. It is a theory for 
dealing with the ground states and their properties'^ . Applications of the Hohenberg- 
Kohn theorem and the Kohn-Sham construction is limited to the ground-state because it 
is the ground-state density of an electronic system that determines the Hamiltonian, and 
consequently other physical observables of the electronic system. This suggests both the 
ground as well as excited-state properties can be determined from the ground-state density 
through the Hamiltonian operator since it characterizes all the states of a system. On the 
other hand, to develop an excited-state DFT akin to it's ground-state counterpart, it is 
important to describe an excited-state of a system in terms of the density of that state. 
Almost for the last two and half decade researchers have investigatect^ i^i^i^'^i^'^i^^i^^i^^i'^^i'^^i'^'''i'^^ 
the possibility of giving a formal foundation to the time-independent excited-state DFT 
like the HKS DFT does for the ground-states, but it is yet to come to its full fruition. 



To make the excited-state calculations feasible within the density functional formalism 
there are two open questions to be answered: (i) Does there exist a mapping between the 
excited-state density and the external potential like the ground-states ? (ii) Secondly, for the 
determination of the excited- state energies, is it possible to construct reasonably accurate 
exchange-correlation energy functionals? We note that although the exact form of the 
exchange-correlation functional is unknown for the ground-states, there are several accurate 
and approximate functionalsi^iSSiSLSSiS^in traditional DFT. The issue is to find such approx- 
imate functionals for the excited-states. The second question is partly answered through an 
attempt made by the present author^ in the recent past by developing an exchange energy 
functional within the local-density approximation(LDA) for a particular class of excited- 
states. This suggests that a correlation functional can also be developed in a similar fashion. 
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We now address the first question which is the main focus of the present study. Over the 
past few years, a lot of attention has been paid to the question of mappin g^^i^^i^^i^^i^^i^° from 
an excited-state density p(r) to the external potential Vextii^, because the entire structure 
of time- independent excited-state DFT depends on that. A brief account of this is as follows : 

The first step in establishing a mapping from an excited-state density to a many 
electron wavefunction is taken by looking for p-stationary statea^^. These are states \E' that 
reproduce a given density (density of the k^'^- excited-state) and simultaneously make the 
expectation value ^\E'|T-|- Kel^^ stationary. However, for a given density there are many 
p-stationary states and thus establishing a mapping requires an additional input. Levy and 
NagyiS*2& provide this by requiring that \E' be orthogonal to ^'^(j < k), which are to be 
determined by the ground-state density po for the system under study. Thus the knowledge 
of the ground-state density is essential to define a map pk — > "^k- This is reasonable 
because it is the ground-state density that really fixes the Hamiltonian of a system uniquely. 
Levy and Nagy then go on to construct a Kohn-Sham system of non-interacting electrons 
whose m^^ excited-state produces the given excited-state density. Further they put forth 
a criterion that the ground-state density of the Kohn-Sham system is closest to the true 
ground state density of the system in the least square sense. 

We have been studying the mapping from an excited-state density to the potential and 
have been investigating the Levy-Nagy criterion for constructing the Kohn-Sham system for 
excited-states. Our work in this direction forms the content of this paper. We show: 

(i) an explicit construction of the external potential from an excited-state density using 
Gorling 's approach^^, 

(ii) that the Levy-Nagy criterio n^^i^^ of constructing the Kohn-Sham system is not perfect. 
We give reasons for it, 

(iii) that for a given state the Kohn-Sham system should be constructed by comparison of 
the kinetic energies of the true and the non-interacting systems. 
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II. GORLING/LEVY-NAGY FORMALISM 



A. Gorling's Stationary-State Formulation 

We start with Gorling's formulation^^ of the excited-state problem. The formulation is 
based on the constrained-search approach^ and provides a map from an excited-state density 
to a many-body wavefunction. For a given excited-state density p{r) a p— stationary state 
\E'(r) is given by making the expectation value ^\I/|T + V^e|^^ stationary with the constraint 
that the many-particle wavefunction \E'(f) giving the same density p(f) . Corresponding to 
each p— stationary state \E'(f^ there is an external potential Vext{^ ■ This has been shown 
by Gorling. We give a different proof of it here, with the external potential Vext{r) arising 
as the Lagrange multiplier to ensure the constraint of producing the density p(f). 

A p— stationary wavefunction \E' gives, by the stationarity principle 



6^\T+Vee\^} + {^\T+Veem}=0 (l) 

with the constraint that 

5p(r) = j {^6^^ + ^^5*} df2 dfN = (2) 

Writting Eq. 1 in expanded form, we get 

j 1^^^ (f + fee) ^ + (r + Vee) } dndf2 rffv = (3) 

Because of condition (Eq. 2) above this will be satisfied if 

[f + V;.) f„),f. ^ / Mr-mr-. r. ... r„)*, . (4) 

Thus in general p— stationarity of \E' implies that it satisfies 

(r + Ke) ^(r-'i,r2,...,r-W) = (J2f^ir^)j ^(fi, r^, rW) . (5) 

However, since T,Vee are symmetric operators and \E'(fi,r2, ...,f/v) is antisymmetric, it's 
necessary (^. fi{ri)) must also be symmetric. Thus all /i's must be the same function /(r). 



Identifying this function as /(r) = —v{r) + E , where hm f (r) = 0, we get \1/ satisfying 

We note that for different p— stationary states Vexti'^ will be different. Thus by applying 
the constrained search method we get many p— stationary states and the corresponding 
external potentials. The question is which one of these corresponds to a given system. 
Levy and Nagy identifyi^*^^ this system as the one where is orthogonal to "^jij < k) 
for a given ground-state density po- Thus in the Levy-Nagy theory''^^, the wavefunction 
\E'[p;po] is a bi-functional of p and po. One subtle point about the Levy-Nagy theory is 
that if the search for is restricted to the space orthogonal to "^jiJ < k), the variational 
principle becomes a minimum principle. The prescription above also makes the functional 
F[p; Po] = min (\E'|r + Kel^) a bi-functional of the excited-state density p(r) as well as 

>I'^p(r) 

the ground-state density po{r) . 




r (a.u.) 

FIG. 1: Shown in the figure are the KS potentials for the excited-state density of ls^2s^2p'^ {^S) 
state of Be generated for the original and one alternative configuration. 

Next we discuss how a Kohn-Sham(KS) system can be realized for an excited-state den- 
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sity. To get a KS system, the mapping from a given excited-state density to a non-interacting 
wavefunction is estabhshed by making the expectation value ^$|T|$^ stationary with re- 
spect to variations in $ with $ giving the excited-state density. How this is done is described 
later. As is the case for the interacting systems, there may be several $'s and the corre- 
sponding KS potentials vks that give rise to an excited-state density and make (^<I>|T|<I>^ 
stationary. Two such potentials for the density of 15^25*^2^^ i^S) of Be are shown in Fig. 1 . 
The excited-state density used is obtained by solving the excited-state KS equation with the 
Harbola-Sahni(HS)^ exchange-only potential. The question that arises again is how do we 
choose one particular KS system to represent a system in its excited-state. An intuitive way 
would be by comparing the ground-state densities, as was done for the interacting systems. 
However, the ground-state density of a non-interacting system that reproduces an excited- 
state density may not be the same as that of the true system (it will be the same only if 
the electron-electron interaction were absent). Thus Levy and Nagy propose^^^ that of the 
many Fermionic non-interacting systems that give the same excited-state density, the one 
whose ground-state density resembles the exact one in the least square sense be identified 
as the Kohn-Sham system. The criterion is obviously exact^° for systems with no electron- 
electron interaction, as stated above. For interacting-electron systems, the criterion appears 
to give'^^ the true Kohn-Sham system, consistent with the orbitals to which the electrons 
have been excited, as we discuss below. 

B. Levy-Nagy Criterion 

We have been investigating the Levy-Nagy criterion for many different excited-states of 
atomic systems. We find that whereas the criterion is satisfied for a large number of excited- 
states, there are exceptions also. We present and discuss these results below. The results 
are obtained as follows. Working within the central field approximation, we perform the 
excited-state calculations using the exchange-only HS method and obtain the excited-state 
density (KLI-OPM^^*^^^ gives similar results). We then use the Zhao- Morrison-Parr (ZMP) 
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method^^ to generate the same density with different configurations and test the Levy-Nagy 
criterion for these configurations. In the ZMP method the Kohn-Sham system for a given 
density p{r) and a configuration of choice^^ is obtained by solving the equation 

^^.ri,|0,(f/)|2-p(f/) 



■dff > 0i(r) = ei(t)i{r) (7) 



2 J \f — f/\ 

in the hmit of A ^ oo . Here {%} are the occupation numbers of the orbitals according to 
the configuration chosen. We choose A large enough so that not only do the densities match 
to a high degree of accuracy, the highest occupied eigenvalues Smax of the above equation also 
matches with the original Smax to within 5%; in fact it is better than within 2% in many 
of the cases. For example, for the ls^2s^2p^ excited-state of Be, we have generated the 
same density with three other configurations: ls^2s°2p^, ls^2s^2p^ and ls^2s^2p^. The Emax 
and the expectation values (-^) , {R) and {R"^) for different configurations are compared 
in Table I. To judge the numerical accuracy of our ZMP program, we also generate the 
excited-state density with the original configuration and compare numbers obtained with 
the original numbers. We see that with the original configuration, the Smax comes to within 
2% of the original value with A = 5000 whereas the various expectation values are essentially 
exact. For the three alternative configurations, the accuracy of , {R) and (R^) is about 
the same but the Emax values differ slightly more depending on the configuration. The worst 
case is the ls^2s^2p^ for which Emax = —0.626 Ry for A = 30000. To make sure that the 
eigenvalue will eventually converge to Emax = —0.658 Ry, we performed calculations for 
different values of A for this configuration and found that Emax = —0.621 Ry, —0.624 Ry 
and —0.626 Ry for A = 10000, 14000 and 30000, respectively, thereby shifting towards the 
true eigenvalue albeit very slowly. We also mention that for the configuration ls'^2s^2p^, the 
uppermost orbital is 2s and not 2p. The local potential in which the electrons are move is 
then given as: 

VKsir) = A / ^ — — Uf/ . (8) 



|r — r/\ 

Having found the potential above, we obtain the ground-state density po{r) of this potential 
by occupying the lowest energy orbitals with the given number of electrons. We calculate 



its mean square distance from the true ground-state density poir) as^ 

A[po(r),po(r)]= / {poW - Po(r)}MV . (9) 

J oo 

By the true ground-state density here, we mean the ground-state density obtained with the 
Harbola-Sahni exchange potentiaL The results for A for a number of atoms and their excited- 
states are shown in Tables II & III . The original configuration of the excited-sates is shown 
in the second column. The third column shows the alternative configurations using which we 
obtain the same density and the fourth column the corresponding A [po{r), po(r)]- It is seen 
from the results that for most of the cases A is the smallest for the original configuration 
but there are cases where A is smallest for a different configuration. For example, there is 
the excited-state ls^2s^2p^ of Be for which the configuration ls^2s^2p^ gives the smallest 
A. Similarly for the state 1522522^3352 ls^2s'^2p^3s'^ of A^e and ls^2s°2p^ of B+, 

ls^2s^2p^ of A^e^"^ A is the smallest for a configuration other than the original configuration 
of the system. Thus we find that the Levy-Nagy criterion, as quantified by Eq. 8 above, is 
not satisfactory in that it leads to a KS system where an excited-state configuration is not 
consistent with the original system. 

C. An Alternative Criterion 

Having found that the Levy-Nagy criterion is not fully satisfactory to identify a KS sys- 
tem, we have looked for other ways of doing so, remaining within the Levy-Nagy proposal 
of comparing the ground-state densities. Thus instead of comparing densities directly, we 
compare them energetically as follows. After obtaining many different non-interacting sys- 
tems for an excited-state density, we take their ground-states and calculate the expectation 
value of the true ground-state Kohn-Sham Hamiltonian(constructed using the HS exchange 
potential) with respect to these ground-states. Thus the calculation proceeds as follows: we 
solve the HS equation for the ground-state of a system and obtain the ground-state Kohn- 
Sham Hamiltonian Hq. The expectation value (Hq) with respect to the true ground-state 
orbitals is designated as {Ho)true'^ the sum of the eigenvalues of the ground-state orbitals. 
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Next we take the different non-interacting systems giving an excited-state density, consider 
their ground-states and calculate the expectation value (i?o)a/i with respect to these states. 
Because of the variational principle (-f^o)a/t should always be above {Hq)^^^^. We then iden- 
tify the true KS system as that for which {Ho)^i^ is closest to {Ho)^^^^. This comparison is 
made in Tables IV & V for the same systems as in Tables II and III with HS density As 
is clear from the table the alternate criterion is better in that the correspondence between 
the original system and the Kohn-Sham system is restored for F,Ne and Ne^'^. However, 
new inconsistencies arise in Al,Si and P'^ although in these cases the difference in the num- 
bers for the original and the alternative configuration is very small. On the other hand, 
the inconsistency in Be and 5+ remains. We note that this criterion is very sensitive to 
the exchange potential. If calculations are done with KLI-OPM exchange potential, the 
inconsistency remains only in Be and B~^ systems. It is clear from the discussion in the two 
sections above that a criterion based on comparison of ground-state densities of excited-state 
Kohn-Sham systems cannot be satisfactory. 

III. PRESENT THEORY 

Given the background above, we now present a consistent theory of excited-states within 
the rubric of density-functional approach. The principal tenets of the theory are: 

(i) There is a straightforward way of mapping an excited-state density p{r) to the cor- 
responding many-electron wavefunction ^(r) or the external potential v^xti^ using the 
p-stationary wavef unctions. The wavefunction depends upon the ground-state density po 
implicitly. 

(ii) The Kohn-Sham system is defined through a comparison of the kinetic energy for 
the excited-states. This avoids any comparison of the ground-state densities which, as seen 
above, doesn't give a satisfactory way of defining a KS system. We now discuss these two 
points one by one. 

To describe the mapping from an excited-state density Pfc(r) to a many-body wavefunc- 
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tion, we take recourse to the constrained-search approach. This gives, as discussed earher, 
many different wavef unctions \I/fc(r) and the corresponding external potential v^^^{r). If in 
addition to the excited-state density we also know the ground-state density po then Vexti.^ 
is uniquely determined by the Hohenberg-Kohni theorem. Thus with the knowledge of po 
it is quite straightforward to select a particular \E' that belongs to a [p, po] combination by 
comparing v^^^{r) with Vexti'f^- Alternatively, we can think of it as finding \E' variationally 
for a [p, Vext] combination because the knowledge of po and Vext is equivalent. Through the 
constrained search above a functional 

F[p;po] = (v|/[p;po]|r + V;e|^[p;Po]) (10) 

is also defined. The prescription above is similar to that of Levy and Nagy2£. but avoids the 
orthogonality condition imposed by them. 



The densities for different excited state for a given ground-state density po or external 
potential Vext can thus be found as follows: take a density and search for \E' that makes 
^|r-|-fee|^^ stationary; check whether the corresponding Vext matches with the given 
Po(or Vext)', if not, take another density and repeat the procedure until the correct p is 
found. Also because of the proof given in section II above, the Euler equation for the 
excited-state density is 

SF[p,po] 

Sp^ff^ +^ext(r1=p (11) 

where p is the Lagrange multiplier to ensure that pfc(r) integrates to the proper number of 
electrons. 

The prescription above for the excited-states in terms of their densities is quite straight- 
forward, particularly because it's development is parallel to that for the ground-states. On 
the other hand, to construct a Kohn-Shan>^ system for a given density is non-trivial; and to 
carry out accurate calculations for excited- states it is of prime importance to construct a 
KS system. Further, a KS system will be meaningful if the orbitals involved in an excitation 
match with the corresponding excitations in the true system. We have shown above that 
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the Kohn-Sham system constructed using the Levy-Nagy criterion fails in this regard. 

In principle, obtaining a Kohn-Sham system is quite easy. Define the non-interacting 
kinetic energy Tg [p, po] and use it to further define the exchange- correlation functional as 

E,, [p, Po] = F[p,po]-lJ J yz^dfdrf - T, [p, po] • (12) 

Then the Euler equation for the excited-state densities will read 

STs[p,po] f pjff) 5E^c[p,Po] , .ION 

+ / -rz ^ar/ H ^— r=r h Vext = A* • (13) 



5p(f) J \f— f/\ ^p(r) 

which is equivalent to solving 

1 



with 



VKs{r) = Vext{r) + / -p — ^rdr/H —r^ — . (15) 

J \r — rf\ op[r) 

However, it is defining Tg [p, po] that is not easy in the excited-state problem. For the 
ground-states, Tg [po] is easily defined as the minimum kinetic energy for a given density 
obtained by occupying the lowest energy orbitals for a non-interacting system. On the 
other hand, for the excited-states it is not clear which orbitals to occupy for a given 
density, particularly because a density can be generated by many different non-interacting 
systems. Levy-Nagy select one of these systems by comparing the ground-state density of 
the excited-state non-interacting system with the true ground-state density. However, this 
criterion is not satisfactory as discussed earlier. Therefore some other criterion has to be 
evolved to construct the excited-state Kohn-Sham system. 



Before searching for other ways of constructing a Kohn-Sham system, we look for 
reasons that may be responsible for the Levy-Nagy criterion not being fully satisfactory. 
We argue that we are not being consistent while comparing the ground-state density 
of an excited-state KS system with the true ground-state density. This is because the 
ground-state density of the excited-state KS system is not the self-consistent ground-state 
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density of the Vextif) obtained for the excited-state density but of a potential different from 
Wext(r) • In comparing the ground-state densities, we are thus not comparing the Vexti"^ of 
the excited-state KS system with the true Vexti'f^, and this at times leads to inconsistent 
results. 



In light of the above remarks, it is important that in constructing the Kohn-Sham 
system only the self-consistently determined quantities corresponding to a given excited- 
state density be compared. Thus we propose that the KS system be so chosen that 
it is energetically very close to the original system. To ensure this, we define the KS 
system as that system for which the non-interacting kinetic energy , obtained through 
constrained search over the non-interacting wavefunctions, is closest to T[p,pq\, already 
obtained through the full constrained search. This then gives the functional Ts[p,pq\. 
Thus Tg [p, pq] is defined as the kinetic energy that is closest to the true kinetic energy 
T [p, Pq] obtained for a given excited-state density p. Defining Tg [p, po] in this manner also 
keeps the DFT exchange-correlation energy close to the conventional quantum mechanical 
exchange-correlation energy. An added advantage of keeping the difference between the 
two kinetic energies Tg [p, po] and T [p, po] smallest is that the structure of the Kohn-Sham 
potential is simple; it is known that contribution of T — Tg gives more structure to the KS 
potential. 



All the statements in this paragraph are justified on the basis of the differential virial 
theorent^i. Using this theorem the exchange-correlation potential for a given density and 
the corresponding many body wavefunction \[' can be written as 



ZKsir; [Fi^^]) - Z(f; [T,]) + J [Vu{f,f/)] [p{f^p{f/) - 2r2(f,f/)] df/ 

-VMr-') = -^ (16) 

p(r) 

where ri(r2) is the first (second) order density matrix and the vector fields Z, Zks are related 
to the true and Kohn-Sham kinetic energy density tensors respectively. For a given [p, po] , 
Z^Ti^T2 are fixed and different configurations chosen give different Z^s- These different 
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configurations can be thought of as arising from a different external potential (as is shown 
below) or from a different exchange-correlation potential^^^. In any case kinetic energy 
difference between the true and Kohn-Sham system isM 



It is this difference that we propose be kept the smallest for the true KS system, and as we 
show below, it gives the Kohn-Sham system consistent with the original system. 

A. Examples 

We now demonstrate the ideas presented above with examples. Since we do not know 
how to perform the general constrained search p , S (^^\T + Vee]"^^ = 0, we take an 
indirect path for the purpose of demonstration. In the following we work with atomic 
excited-state densities generated by HS exchange-only potential and take these as the 
excited-state density. The densities and the energies obtained by the HS formalism are 
essentially the same as those of Hartree-Fock(HF)^ theory for both the ground as well 
as the excited-states. Similarly the HS exchange potential is very close to the true local 
exchange potential of the optimized potential method(OPM). Thus the formalism is well 
suited to test the ideas presented above. Thus we start with this given [p, Vext] combination. 

We first demonstrate that an excited-state density is generated by different external 
potentials depending on the configurations chosen to generate the density. For this we use 
the ZMP method to get the non-interacting potential giving the density p{r) and subtract 
from it the Coulomb fcoui(r) ^^'^ exchange v^^{r) potential to get the external potential 
Vextii^ (^-e. Vextii^ = VKsi^^ ~ "f^couiif) — ^x^if)) ■ The cxchauge potential for a given set of 
occupied orbitals is obtained using the Harbola-Sahni formula^ for it. Shown in Fig. 2 are 
different external potentials thus generated for the ls^2s^2p^ i^S) of density Be. To check 
our consistency, we first obtain Vext{'f^ for the original configuration (ls^2s°2p^) and find it 
correctly to be — - Ry . The other configurations that we use to obtain the same density 




(17) 
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FIG. 2: Shown in the figure are the external potentials for the excited-state density of 
\s^2s^2p^ (^S) state of Be corresponding to the original and one alternative configuration. 

are ls^2s°2p^, ls^2s^2p^, ls^2s^2p^. We have shown only two potentials corresponding to 
the configurations ls^2s^2p^ and ls^2s°2p^, and compared them with the true external 
potential Vext — ~^ ^V- As discussed earlier, only one of these - that corresponding to the 
original configuration - matches with the true external potential. 

Next we show that the external potentials corresponding to the ground-state density 
of excited-state KS systems are different from the external potential for the excited-state 
density. This, as pointed out earlier, sometimes leads to non-satisfaction of the Levy-Nagy 
criterion. Shown in Fig. 3 are the Vexti'f^ corresponding to the ground-state densities of 
different configurations for the excited-state density of ls^2s'^2p'^ {^^S) state of Be. These 
potential are also obtained by subtracting from the Kohn-Sham potential vks the Coulomb 
and the Harbola-Sahni exchange potential calculated by occupying the corresponding KS 
orbitals in the ground-state configuration. We see that whereas the true external potential 
is — ^ Ry, the external potentials corresponding to the ground-states are different. It is 
this difference that leads to the ground-state densities different from the true one, and also 
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FIG. 3: Shown in the figure are the external potentials corresponding to the ground-state densities 
of excited-state KS systems for ls^2s^2p^ (^S) state of Be. The potentials are compared with the 



sometimes to inconsistencies between the KS and the true systems. 

Finally we show that a comparison of the kinetic-energies(KE) leads to the appropriate 
Kohn-Sham system. We demonstrate this with the densities generated using the HS 
exchange potential and with one example with correlated density. The numbers for the 
noninteracting kinetic energy for different configurations corresponding to the same excited- 
state densities as considered earlier are shown in Tables VI and VII and are compared with 
their original kinetic energy. Since the HS potential is a local potential itself, the correct 
configuration gives AT^ = (slight difference arising due to numerical calculations) and 
wrong ones a larger value, as is evident from the numbers in the Tables. We see that unlike 
the Levy-Nagy criterion, a comparison of excited-state KE leads to a proper KS system in 
all the cases. Of course when we use the correlated densities, the difference AT is not going 
to be zero for the proper configuration but should be the smallest. This is clearly seen in 
the example of ls2s(^S) state of He atom where the true T[p] is 2.146 a.u.. We have used 



true external potential Vgxt = ~f ■ 
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the ZMP procedure to obtain the KS potentials in this case also and see that KE for the 
Is^ and the ls2s configurations is 2.044 a.u. and 2.153 a.u., respectively. In the latter case 
it is closer to the true kinetic energy of the system. Thus the configuration ls2s represents 
the KS system for the 1525(^5*) He density. 



IV. CONCLUDING REMARKS 

Building on the work of Gorling and Levy and Nagy, we have presented in this paper 
a consistent theory of excited states within the density-functional formalism. The theory 
is based on constrained-search and defines a bi-density functional F[p, po] without the 
orthogonality constraint of Levy-Nagy. Further the theory gives a clear definition of 
the excited-state Kohn-Sham systems as that whose kinetic and exchange-correlation 
energy components are closest to those of the true system. This avoids the problem of 
comparing the non-self-consistent ground-state densities, as proposed in the LN theory, so 
no inconsistency arises in identifying an excited-state Kohn-Sham system. 

To conclude, we have analyzed the theoretical foundations of excited-state time- 
independent density-functional theory and have put Kohn-Sham calculations for excited- 
states on a firm footing. It is clear from our present work that an excited-state configura- 
tion of the KS system corresponds to a similar excited-state of the true system (with the 
major component of the excited-state wavefunction involving the same orbitals as the KS 
system). The work should also help in providing guidance in construction of the excited- 
state exchange-correlation energy functionals to facilitate self-consistent determination of 
the excited-state energies. 
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TABLE L Shown in the table £max and expectation values (■^) , {R) Sz {R"^) for various configura- 
tions giving the same density as that of ls^2s°2p^;^ S state of Be. The self-consistently determined 
values of these physical quantities are : emax = —0.658 Ry,(^^) = 5.818, (i?) = 6.755&(i2^) = 
17.309. 



configurations 


A 


£max{Ry-) 


a) 


{R) 






5000 


-0.649 


5.818 


6.755 


17.312 


Is22s°2p2 


5000 


-0.649 


5.819 


6.755 


17.312 




5000 


-0.655 


5.819 


6.755 


17.312 




30,000 


-0.626 


5.818 


6.755 


17.310 



19 



TABLE II: Value of A(See Eq. 8) for different configurations (S*" column) giving the same excited- 
state density as that for the original configuration(2"'^ column). Systems where the LN criterion 
is not satisfactory are indicated with a '*' against them. 



atoms/ions 


true.config. 


alt.config. 


A[po(r),po(r)] 






ls^2sV 


0.0662 


*Be 




ls^2sV 


0.1627 






Is22si2pi 


0.8758 






Is22s°2p3 


0.0002 


B 




ls'^2s'^2p'^ 


0.0065 






ls'^2s'^2p^ 


0.0286 






ls'^2s^2p^ 


0.0008 


C 




Is22s22p2 


0.0156 






ls22s°2/ 


0.0181 






ls^2s°2/ 


0.2903 


N 


ls^2s°2/ 


ls22s°2/ 


6.8461 






ls'^2s^2p'^ 


8.8409 






ls^2s^2p^ 


0.3875 


O 


ls^2s^2p^ 


ls'^2sV 


8.9609 






l.s-2.s"J-2/y^ 


12.1209 






ls22s^2/ 


0.0002 


F 




ls'^2s'^2p^ 


0.0621 






ls^2s^2p^3s^ 


0.0521 


* p 


ls'^2s'^2p^3s'^ 


ls'^2s'^2p^3s'^ 


0.2704 






ls^2s'^2p^3s^ 


0.7138 


*Ne 


ls^2s^2/3s2 


ls^2sV^s^ 


1.6176 
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TABLE III: Caption is the same as in Table II. 



atoms/ions 


true.config. 


alt.config. 


A[po(r),/9o(r)] 






ls^2s^2p'^ 


0.0402 


*B+ 


ls^2s°2/ 


l5^25°2p3 


0.1931 






15225^2^1 


2.0261 






15^25^2^2 


0.1465 






l5^25°2p3 


0.3463 






Is22s22p63s°3p2 


0.0007 


Mg 




1522522/35^3/ 


0.0014 








0.0039 






Is22522/35°3p3 


0.0008 


Al 


Is22s22/3s°3j93 


Is22s22p63si3p2 


0.0022 






Is22s22/3s23pi 


0.0073 






1522522/35^3/ 


0.0007 


Si 




1522522/3523/ 


0.0027 






1522522/35^3/ 


0.0078 






Is22s22/3s°3p3 


0.0004 


Si+ 




1522522/35^3/ 


0.0042 






1522522/3523/ 


0.0148 






I522522/35I3/ 


0.0005 


P+ 


1522522/35^3^3 


I522s22/3s23p2 


0.0053 






Is22s22/3s03p'^ 


0.0099 






ls22522/35°3/ 


0.0006 


P 


l522522/35°3/ 


1522522/3513/ 


0.0055 






1522522/3523/ 


0.0207 
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TABLE IV: Comparison of the expectation value —{Hojait. and —{Hq) exact for various config- 
urations corresponding to a given excited-state density. The first column mentioned are the 
atoms/ions, the second column the original configuration , fourth column the alternative con- 
figurations considered. In the third and last column given the expectations values — (-ffo)ea;act and 
-{Ho) alt. respectively. 



atoms/ions 


exact config. 


— {Hq) exact a-U. 


alt. config. 


-{Ho) alt. a.u. 


*Be 




8.9272 


Is22si2pi 
l.si2.s°2/ 
ls^2s^2p^ 


7.4302 
8.8051 
8.8601 


Be 




8.9237 


ls^2s^2p2 
ls'^2s^2p^ 


6.1351 
8.9217 


B 




15.1482 


ls^2s^2p^ 
ls^2s^2p^ 
ls'^2s^2p'^ 


15.0828 
15.1471 
15.1365 


C 




22.9759 


Is22s22p2 

ls'^2s^2p^ 
Is22s02/ 


22.9556 
22.9723 

22.9473 


N 




32.6951 


ls'^2s'^2p^ 

15^25^2/ 

ls22s°2/ 


32.6682 
32.6949 
32.6712 


O 


ls22s°2/ 


43.3618 


Is22s22p4 

ls22s°2/ 
ls'^2s^2p^ 


43.2215 
43.3608 
43.3252 


F 


ls^2s^2p^3s^ 


55.8686 


ls^2s^2p^3s'^ 
ls''2s^2p^Ss^ 


55.4295 
55.6436 


Ne 


ls^2s^2p^3s'^ 


70.1743 


ls^2s22/3s^ 
Is^2s^2p63s2 


69.2101 
69.4399 
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TABLE V: Caption is the same as in Table IV. 



atoms /ions 


exact config. 


— {Hq) exact a-U. 


alt. config. 


-{Ho) alt. a.u. 


*B+ 




16.6626 


Isi2si2p2 
ls^2s°2p3 

ls'^2s^2p^ 


16.6106 
16.5501 
14.2959 






93.0672 


ls^2s^2p^ 
lsi2s°2/ 


92.9276 
92.9699 


Mg 


Is22s22/3s03^2 


110.9422 


15^2522/3523/ 
I522522/35O3/ 

l.s22.s-2/)'^3.si3/)^ 


110.9091 
110.9368 

110.9341 


*Al 




137.1905 


Is22s22/3si3p2 

Is22.s22/3s03p'^ 
Is22s22/3s23pi 


137.1856 

137.1853 
137.1655 


*Si 




166.2441 


Is22s22p63s23p2 

Is22s22/3s^3/ 
Is22s22p63s03p4 


166.2411 
166.2393 
166.2201 


Si+ 


Is22s22/3s°3p3 


170.5966 


Is22s22/3s23pi 
Is22s22/3s03p3 

ls'^2s'^2p^3s^3p^ 


170.5714 
170.5941 
170.5923 


*p+ 


ls'^2s'^2p^3s^3p^ 


203.2722 


Is22s22p63s23p2 
Is22s22p63si3p3 

Is22s22p63s03/ 


203.2693 
203.2687 

203.2522 


P 


Is22s22/3s°3p5 


198.1074 


Is22s22/3s23p3 
Is22s22/3s°3p^ 
Is22s22p63si3p^ 


198.0818 
198.1058 
198.1025 


s 


Is22s22/3s°3/ 


232.1008 


Is22s22/3.s23/ 
Is22s22p63s03p6 
Is22s22p63si3p5 


232.0761 
232.0996 
232.0957 
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TABLE VI: Shown in the table atoms/ions with the original excited-state configuration in the 
second column and density of this generated by various alternative configurations in the fourth 
column. In the third and fifth column given are the values of the kinetic energies corresponding to 
the original and alternative configurations respectively 



atoms/ions 


true config. 


T\p] a.u. 


alt. config. 


Ts[p\ a.u. 


Be 




10.1489 


ls^2sV 

ls^2s^2p^ 
ls'^2s^2p^ 


10.0177 

10.1481 
8.1357 


B 




24.1249 


ls'^2s'^2p^ 
ls'^2sV 
ls'^2sV 


23.7627 
24.1211 
23.9238 


C 




37.5938 


lsHs'^2p'^ 
ls^2s^2p^ 
ls'^2s^2p'^ 


37.2985 
37.5922 
37.9299 


N 




38.5551 


ls'^2s°2p^ 
ls^2s°2p^ 

ls^2s^2p'^ 


30.5856 
38.5525 
30.6238 


O 


ls^2s^2p^ 


54.7136 


ls'^2s^2p^ 
ls^2s^2p^ 
ls'^2sV 


44.6244 
54.7095 
44.5899 


F 




98.5267 


ls'^2s'^2p^ 
ls^2s^2p^ 


97.8733 

98.5212 


F 


ls^2s^2p^3s'' 


98.2631 


ls'^2s'^2p'^3s^ 
Is22s22p33s2 


97.8746 
98.2393 


Ne 


Isi2s^2/3s2 


94.6521 


ls^2s'^2p^3s^ 
ls^2s^2/3s2 


93.4337 

94.6364 
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TABLE VII: Caption is the same as in Table VI. 



atoms/ions 


true con fig. 


T[p] a.u. 


alt. config. 


Ts[p] a.u. 


B+ 




16.8390 


15^25^2^2 

ls^2sV 
ls'^2s^2p^ 


16.6581 
16.8378 
13.2701 






76.5893 


ls^2s^2p'^ 
ls^2s^2p^ 


76.3318 
76.5837 


Mg 




199.3771 


ls'^2s^2p'^3s^3p^ 
Is22s22/3s°3p2 

l.s22,s22/,fi3.s23/>" 


199.2404 
199.3661 

199.1455 


Al 




241.5112 


Is22s22p63si3p2 

ls^2s^2p^'3s°3p^ 
ls^2s^2p^3s^3p^ 


241.3098 

241.4967 
241.1428 


Si 




288.7507 


lsHs'^2p^3s'^3p^ 
ls^2s'^2p^3s^3p^ 
ls'^2s'^2p^3s^3p'^ 


288.4802 
288.7335 
288.9977 


Si+ 


ls^2s^2p^^s^'ip^ 


288.0463 


Is22s22p63si3p2 
ls^2s^2p^3s^3p^ 
ls'^2s'^2p^3sHp^ 


287.7594 
288.0259 
287.5146 


P+ 


Is22s22p63si3p3 


340.1993 


lsHs''2p^3sHp'^ 
Is22s22p63si3p3 

ls'^2.s'^2p^3s^3p'^ 


339.8495 
340.1796 

340.5202 


P 


Is22s22p63s03^5 


339.8574 


lsHs'^2p^3s^3p'^ 
Is22s22/3s°3p5 
lsHs'^2p^3sHp^ 


339.5241 
339.8390 
339.2355 
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